Irreversibility from a Reversible Equation 
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After a discussion on the state of local equilibrium 
with temperature inhomogeneity, comparing mixture state 
reprsentation in statistical mechanics and pure state repre- 
sentation in thermo field dynamics, a simple model is solved 
to show that a reversible equation of motion with the initial 
condition having inhomogeneous temperature can lead to ir- 
reversible, viz. diffusive, behaviour. Yet, the solution is time 
symmetric exhibiting diffusion both towards future and past. 



I. INTRODUCTION 

The irreversible processes remain to be understood de- 
spite the long history and many attempts ||. In 
particular, those processes caused by temperature gradi- 
ent, say, heat conduction and thermal diffusion, have no 
approaches from the first principles that are commonly 
accepted, while those driven mechanically, say, electric 
conduction, has been handled by the perturbation the- 
ory in Hamiltonian dynamics, at least for the purpose of 
computing the transport coefficients ||. 

We have been trying to obtain a microscopic pictures 
for the heat conduction (J] and thermal diffusion || by 
using the thermo field dynamics Q. However, the cal- 
culations got so involved as to obscure the picture. We 
are led to try a simple model, if not quite realistic, to see 
in a mathematically transparent way that an irreversibil- 
ity can result from a reversible equation when the initial 
condition involves temperature inhomogeneity. In this 
calculation, we use the standard machinery of statistical 
mechanics. 

After a brief introduction to the thermo field dynam- 
ics in Sec. 2, we shall relate in Sec. 3 the formulations of 
the state of local equilibrium under temperature gradi- 
ent in the thermo field dynamics and in the more usual 
statistical mechanics. It is interesting that the thermo 
field dynamics exhibits a non-local feature, though micro- 
scopic, of the local equilibrium states; statistical mechan- 
ics gives an equivalent description but in a more indirect 
way. Then in Sec. 4, we employ a simple model to show 
in a mathematically transparent way that a reversible 
equation of motion with the initial condition having in- 



homogeneous temperature can lead to irreversibility. Our 
result is irreversible in that it shows a diffusive behaviour 
among other features, and yet symmetric with respect to 
the direction of time in that the diffusion towards fu- 
ture/past comes out when the equation is solved towards 
future/past, respectively. 



II. THERMO FIELD DYNAMICS 

The basic idea of the thermo field dynamics is to repre- 
sent the thermal state not by the mixture state with the 
density matrix p((3) but by a pure state |/3) which can be 
obtained from p((3) by the GNS construction. Thus, the 
thermal average of an observable A reads 



Tr[p(f3)A} = (l3\A\f3), 



(1) 



the left-hand side belonging to the usual statistical me- 
chanics and the right to the thermo field dynamics. 

With any complete set {\n)} of orthonormal states, one 
can write 

|/3>:=$>(/?) 1/2 |n)®<n|. (2) 

n 

In fact 

m n 

= Y,{n\pA\n) =Tr[p(p)A}. 

n 

The thermo field dynamics further introduces opera- 
tors that act on the bra part as 



(I ® A)\iP) ® ( X \ :=|V>® (^Ix)) 1 



(3) 



thus doubling the number of operators to be considered. 
One uses shorthand notations, 



A*:=I®A, |V) (xl 
Then, one finds: cA — c* A, and 



ip) ® (x|- 



(4) 



1 



[a,a*] = l => [5, 5^ = 1, etc.. (5) 
Suppose a acts on \n) as an annihilation operator, 
a\n) = y/n\n — 1) 
and p = Z^ 1 e~ l3n with H — wa'a, so that 
| / 3> = Z 1 (/3)- 1 /2^ e -^/2| n)(n | j 



(6) 



then 



^ Ee -*/^ ln _ 1)( , 

v 7 n— 1 



-/3lo/2 



a 



1 OO 

___ £ e -»/W2V£ + T|n><n + 1|, 



implying 



(a-e-^ 2 St)| 



= 0, 



which is called the thermal state condition. In terms of 
9 as defined for the f3 by tanh6> = e~ 0w / 2 , 

(a cosh 0- aUmh6\\l3) = 0, (7) 

which reminds us of the Bogoliubov transformation (?]] : 
UiP)aU(P)- 1 = acosh0-S t sinh0, (8) 

with 



U(f3) — exp (^a* a* — aaj 



This operator generates the state \0) from the "no par- 
ticle state" |fi) := |0)(0| of a, 

U(J3)\n) = 1/3) (9) 
as one sees from the operator (||) annhilating (^J) , or from 

U(m) = exp^attanhflp) 

1 OO 

= ^7mT7jE e ""^ /2 NH = l/3>- (10) 

n =0 

We have so far considered only a single mode. In gen- 
eral, however, the object has many modes, and often in- 
finitely many, to which a tensor product of the states 
considered above corresponds. 



III. NONUNIFORM TEMPERATURE 

Suppose we change the temperature a bit uniformly all 
over the system. Then, 



\f3 + A(3) = U(f3 + A(3)\Q) 
1 



exp 



Q t 5 t e -(/3+^W* 



Z 1 {[3 + Apy/ 
which can be approximated as 



\[3 + A(3)=N 1 (x-^e-^/Vst 



l«>. 



with N = Z 1 (/3) 1 ^ 2 /Z 1 (f3 + A/3)^ 2 , the normalization 
factor. 

For a system with many modes labelled by k, Z\ is 
replaced by Z, and 



with 



Define: 



\/3 + A(3)=N(i-ir)\l3) 



— - — e k/ a k 'at '. 

k 



(11) 



^ V k 



0*0 ■= -^/Z ake 



ik-: 



(12) 



then 



G = -iA0 J V (x+ Q <7(O0 f f sb - |) (13) 

For a gas of free particles (or quasi-particles), for which 

k 2 



2m' 



(14) 



one has, in high temperature approximation, 

,3 \ 1/2 



.9(1) 



m 



3 exp 

(3h 2 1 



which is short-ranged as shown in Table 1. 



(15) 
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Since <?(£) is short-ranged, the operator ( |l3| ) can be re- 
garded as a sum of local operators. Therefore, the state 
with nonuniform temperature (3 + f3 a (x) may be repre- 
sented as 



\P + Pa) =N(l-iT[P a ] 



(16) 



with the temperature inhomogeneity /3 a (x) being pushed 
into the integrand of (Jllf): 

T\j3 a ] =-iJ<p(x+ (x - |) g(^)f3 a (x)d 3 xd 3 ^ 

(17) 

since the temperature inhomogeneity (3 a (x) we are con- 
cerned with in studying thermal processes varies only 
slowly on a macroscopic scale. 

We add that ([131) is often written as 



r/ ; , = / <i»i ( x + f\ r 2 $ (x- | ). <jii :ij„u- )•!■>. ni \. 



in terms of the thermal doublet, 



<f>(x) 



i 



(18) 



(19) 



This is the representation we have been using for the 
state with nonuniform temperature in our studies of heat 
conduction and thermal diffusion Q| ||. 

As is clear from its derivation, this representation is 
equivalent to the one by the density matrix, 



p[P + P a {x)} := -exp 



-0H 



(3 a (x)H(x)d 3 x 



(20) 



or, more precisely, to its expansion to the first order in 

A*, 

p[(3 + (3 a {x)] = Z- 1 e- ' H 

x ^1 - J d\e xm J d 3 xl3 a (x)H{x)e- xl3H ^j , (21) 

provided that (3 a (x) is slowly varying. 

It is interesting to compare ( pi] ) and ( |l6| ) + (17) ob- 



serving that the effect of the temperature inhomogeneity 
f3 a (x) appears not to be strictly local in thermo field dy- 
namics but has a small extention as represented by <?(£). 
We repeat that the two repesentations are equivalent, 
and the effect of g(x) is also there in ( f2l| ) as a little anal- 
ysis shows. It is only that the effect can be seen more 
easily in the thermo field dynamics. 

Another remark concerns Zubarev's assertion j^] that 
no irreversible processes can follow from the density ma- 
trix (|20|). Our model calculation in the next section will 
provide a counterexample to his assertion. Further dis- 
cussion will be given in Sec. V. 



IV. A MODEL EXHIBITING IRREVERSIBILITY 

To sec that irreversibility can result from the interac- 
tion of particles even when their equation of motion is 
reversible, we consider a simple model B Q of a Bose 
gas (particle mass: 1/2) with a quadratic Hamiltonian, 



H 



^ VopOp + y X! v o{alaL p + 2aj 3 a p + a p a_ p ). 



p 



p<K 



(22) 



This Hamiltonian is time-reversible, as one sees by apply- 
ing time reversal p — ► — p to the momenta that label the 
operators et p . The price we have to pay for the simplicity 
of its analysis is the non-conservation of the particle num- 
ber; this is serious because we are going to look at the 
particle diffusion as a sign of irreversibility. Nevertheless, 
we shall be able to see this sign |ll[, though superposed 
with an unwanted decay of theparticle number. 

We use the representation ( P0|) of the temperature in- 
homogeneity in terms of the density matrix, for which we 
need the Hamiltonian density, 



H(x) 



d(fv(x) d(f)(x) 
dxi dxi 

+y J v(r) 4>\x + ~r)4t(x - ^r) 

+(j>>(x + ^r)<j)(x - -r) +4> ] {x- ^r)<j)(x + ^r) 



+cj){x + -r)cj)(x - ir) 
where v(r) is such that 



v(k) = / v(r)e 



d 3 i 



d 3 r 



vo 




(k < K) 
(k > K), 



(23) 



v(r) = A^voK 2 ji(Kr) jr. 
In terms of the Fourier transform, we find 



v(r)${x + -r)(f)(x — —r)d 3 i 



= ^E-(^)4'«p^ p ' ) ' x 



and similarly for their Hermitian relatives. We assume 
for simplicity that K is sufficiently large to make all the 
v{k)'s for thermally agitated particles to have the con- 
stant value vq; we shall hereafter suppress the proviso 
P < K . 

For the particle number density, 

io(x) = tf(x)<f>(x) = I Y, «p'ape l(p " p,) - x , (24) 
p,p' 
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we wish to find out the long-time behaviour of the devi- 
ation, 



(jo(M)> {1) 
1 



Tr 



m I dX d 3 x'[3 a {x')H{x') 



from the average, 



Oo(.x,t))^ = ^Tt[e-^j (x)] 



(26) 



where 

t = t + i\/3 
and we have assumed that 

Pa(x)d 3 X = 0. 

The Hamiltonian ( |22| ) can be diagonalized as 

H = ^ w p q p«p ( w p : = P\/p 2 + 2c 2 , c 2 := po^o) , 
p 

(27) 

where the zero-point energy has been dropped, by a Bo- 
goliubov transformation 

t 



a p = a p cosh 8 p — ac p sinh p , 



with 



sinh 8 V 



p 2 + c 2 



2p v / p 2 + 2c 2 2' 



(28) 



cosh 6„ sinh 8„ = , . 

P P 2 P yJp 2 + 2c 2 



The thermal average (25) turns out to be given by the 



(.7oOM)> (1) = (joOM)> (8 ' _) + (j (x,t))^ 

+(j ( x ,t))^ + Oo(*,*)> (c,+) , 



where 



(joOM)> (k,±) 
-1 



2(2tt) 6 



d 3 k/3 a (fe)e lkx / F^flfe.A) (29) 



with 



F( s '±)(fc,r) = f d 3 Pg(p',p)smh(9 p ,+6 p ) 

f /(a; P 0/K)e- i ^' + ^) r r , m 
X \{/K) + l}(/K) + l)e^' + ^ ( ' Uj 



and 



F< c,± >(fc,T) = /" d 3 Fh{p',p) cosh(<V 



r /Ko{/K) + i}e-^-^)- 

\ {/K0 + i}/K)e 4( ^'-^ )r , 

(25) Here, p' = P + fe/2, p = P fc/2, and 



(31) 



g(p',p) 
h{p', P ) 



(p'-p) 



sinh 
cosh 



=Fc (cosh^p/ — sinh p / ) (cosh P — sinh6>p), 



(32) 



of which g is associated with the "off-diagonal" term 
a p' a -p an< ^ ^ w ^h the "diagonal" a p ,a p in the Hamil- 
tonian density TL(x) . We note 



Tr 



= /(Wp')/(^p)((5p' _q'<5-p, q + <5p',q<5p,q')- 



(33) 



A. Evaluation of the integral f( s ± ' 

As a key for the following calculations, the tempera- 
ture inhomogeneity /3 a (x) is assumed to be varying on 
a macroscopic scale, so that the support of its Fourier 
transform is very small in comparison with the scales on 
which ojp and the coefficients of the Bogoliubov transfor- 
mation vary. 

We begin the evaluation of ( p9| ) from the P- 
integration, 



F (s ' + \k, T )= / P 2 dP sinh[20 F ] (7(w F ) 



x / dn P g(p',p)e-^'+^\ (34) 



where, in sinh and /, we have put p = p' = P taking 
advantage of the small support of (3(k). However, we can- 
not do so in g and the exponent, because g would vanish 
if we did and the exponent is proportional to r that in- 
volves the macroscopically large variable t . Instead, we 
expand to order k 2 , 



tjJ v i + LU V 



2lu t , 



P 2 



2LJ r : 



-(P.fc) 



(35) 



and put 



9(p',p)=9o(P)+9i(P)(P-k) 

<P ■ k) 2 

+g 2 (P)k-k + g 3 (P) [ -^- 

finding, however, that go — g(P,P) = by the Bo- 
goliubov transformation diagonalizing Tl [cf. the remark 
after @], g x = by the symmetry g(p',p) = g(p,p') 
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from (|32j) ; we take in advance the later result into account 
that 

g 3 (P) = cosh26> P sinh 2 29 P 
vanishes at the value of P that concerns us. Thus, 

9(P',P)=92(P) with 52 (P) = -isinh20 P . (36) 

The integration over the angles in (|34|), 

2 



Io= I < 
/o 



MCOS 7 sin 7 d 7 = - / dy, 
<* Jo 



is carried out by adding and subtracting 



e y dy 



taking advantage of Im a = A/3 > 0, where 7 is the angle 
between P and k, 



c 4 P< 



k 2 r, a := (— ia) 1 ^ 2 , and y — ax. (37) 



We note — 7r/2 < arg a < and |ct| 
Then, asymptotically 



00 as t 



Thus, 



V* 1 



F (s ' +) (fc,T) 

-2vr / P 2 dP sinh[26»p](/(£jp)) 3 2 (P)fc 



(38) 



V^F 1 



■ exp 



,c 4 P 2 



•fc 2 T 



x exp 



i [ 2uip + 



P 2 



2cjp 



(39) 



Similar calculations give 



2ir I P 2 dP smh[26p]^f(oj P ) + 1 ) g 2 (P)k 2 
.c 4 P 2 



A 1 
x <; 1 h exp 



-k 2 i 



x exp 



i 2wp 



P 2 



2cop 



-k" 



(40) 



If we change the variable of integration P to — P in j4C| ) 
after analytically continuing the integrand to the complex 
P-plane with two cuts to cope with the singularity of 
(P 2 + 2c 2 ) 1 / 2 ; one extending along the imaginary axis 
from i\[2c to ioo, and the other from —iy2c to — ioa, 
then cj_p = —tup, and consequently 



f{u- P ) + 1 



e -(3Ku P _ ^ 

sinh0_p = — sinh^p 



= -f(cop), 



(41) 



so that for P< s ) := p( s >+) + 



I" dP G(P) ( ^ - -1 exp [i^fcVl \ 

J-OO I O- fJ 2 L J J 



x exp 



P 2 + c 2 

-i ( 2wp H A; 2 I ; 

2uop 



(42) 



where 



G(P) := 2^P 2 sinh[20 P ]{/( W p)} 2 ff2 (P) k 2 



_Z. Long-time behavior of 

The time scale we are interested in is macroscopic. To 
evaluate the integral (|42| ) asymptotically for t — ► 00, we 
use the saddle-point method. Candidates for the saddle 
point P s .i for the part of the integral 



Ji := 



x exp 



1 



i < 2ujp 



dPG(P) 

P 2 + c 2 
2lo p 



c 4 P 2 



(43) 



are determined by the condition of vanishing derivative 
of the exponent. 



dlUf 



dP 



2u 2 p 



c 4p2x, 


duip 


LUp J 


dP 


' P 


2c 4 P 


Wp 


LUp 



k 2 = 0, 



where 



dLUp 2{P 2 



dP V-P 2 + 2c 2 

For small fc, the main part of the root is given by P = 
±ic. Then, we put P 2 = -c 2 + £, finding £ = -(3/4)/c 2 
up to order fc 2 , and hence P S)1 = ±ic{l + (3/8)(£; 2 /c 2 )} 
as candidates for the saddle points. Since 

d 2 uu P 2P(P 2 + 3c 2 ) 



dP 2 (P 2 + 2c 2 ) 3 / 2 

takes the value ±4« + 0(k 2 /c 2 ) at the saddle points, the 
exponent behaves like 



2uu f 



c 4 P' 



2tu f 



k 2 



= ±i{(2c 2 + k 2 ) + 4(P - P M ) 2 } 
in the neighborhod of P s ,i, so that 

J, ~ _ J_ e ±( 2 C 2 +fe 2 )- / G(P)e ±4 ( p - p - 1 ) 2r dP (44) 
ct 2 J-00 
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For r ~ t — > +00, therefore, we should take the lower 
sign here 



1 k 2 



(45) 



for the saddle point which the path of integration should 
pass. Thus, 

F(S) As2G( p s,i)e- {2c2+k2)T / e" 4 " 2 d s , (46) 

and we have only to carry out the expansion ( |36"1) around 
P = P s i , obtaining 

^«(fc 1 r) 1 ~-i(I) 8/a /(a, ft . 1 )V^+^. (47) 

Note that the small factor k 2 of G has been cancelled by 
1/ct 2 . The factor e~ k T is a welcome sign of the diffu- 
sive behavior. However, it is accompanied by e 
damping, which is an unexpected feature. 



-2c" T 



,2. Long-time behavior of F 2 
We still have to evaluate the other part of the integral 



x exp 



dPG(P) 



i [ 2u>p + 



P 2 + c 2 
2wp 



(48) 



Candidates for the saddle point P S; 2 for this integral are 
given by the roots of 

'P 2 + c 2 duj P P 



duip 



and hence 



2lo 2 p dP lup 



1 k' 



k 2 = 0, 



P s , 2 =±ic[ 1 + - — 



±ic 1 



2\ V2 



up to order k 2 . Since 

8c 2 y v 4 

we see the exponent in (ktSl) behaves as 



(49) 



±ic + Q(k 4 ), 



i [ 2ujp + 



P 2 



2u!p 



k 2 )T = ±2{c 2 + 2{P-P s . 2 ) 2 }T 



in a neighbood of P s ,2- Therefore, we choose the saddle 
point P Si 2 = —ic I 1 + — — 5- I to let the path of integra- 



2 e- 2c \ 



tion ( |48| ) pass, obtaining 

which implies a simple damping without any diffusive 
character. 



(50) 



3. Integration over X 



It remains to integrate (47) and (|50| ) over A. We may 
consider only the part, 



-He' 



^d\ 



1 /3c 



cot- (51) 



2/3c 2 2 



Thus, summarizing 



3/2 



-2cHt / / -Art 



V^kH) f3 a (k)d\. (52) 



Two remarks are in order. (1) e~ k 1 leads to a dif- 
fusion. This does not contradict the basic reversibility 
of our model. If we reversed the time to look at the 
evolution towards the past, t < 0, then we must take 
the path of integration crossing the other saddle point as 
given by (|4j]) with the upper sign, obtaining the diffusion 
towards the past. Our result is time symmetric in this 
sense, yet exhibiting the irreversible features. The same 
circumstance (but no damping) was found in our previ- 
ous approach jq] with the thermo field dynamics using a 
more realistic model. 

(2) The simple damping implied by ( |50| ) should call for 
an interpretation if it could be taken as an irreversible 
behavior. Its amplitude has a factor k\/t extra to Ji7|), 
k being small but t large. The damping might be an 
artifact due to our model Hamiltonian not conserving 
the particle number. 



B. Evaluation of the integral (31) 

Since 

h{P,P)=u> Pl cosh[20 P ] 
P 2 + c 2 



P 2 + c 2 



w P+k/2 — w P+k/2 



LOp 



UJp 

-(P-k) 



(53) 
(54) 



( |3l| ) turns out to be 

poo 

p(c±)= / P 2 dP(P 2 +c 2 )f(ojp){f(oj P ) + l} 
Jo 

r r p 2 + c 2 

X / dQp exp T?? ■ fcrrnsTi . 

J 1_ VP 2 + 2c 2 . 

The integration over angles gives the same result /2 
for both p( c '±) : 

ip(c) 
2 



,/ P 2 dP^l 

ikr 



f(u P ){f(u;p) + i} 



J^zbexp 



p2 , 2 

±2i + : fcr 
VP 2 + 2c 2 



(55) 



G 



We use the saddle-point method after analytically con- 
tinuing the integrand on the complex P-plane with two 
cuts, [iV2c, oo) and (— oo, — i\/2c] as before. The saddle 
points are determined from the exponent by 



P 2 



2 _ P(P 2 + 3c 2 ) 

dP {P 2 + 2c 2 ) 1 /2 _ (p2 + 2c 2 ) 3 / 2 







(56) 



as P s = ±i\/3c + e, two twins, each of which, being sep- 
arated by one of the cuts, consists of the right and the 
left member with an infinitesimal e = +0 and e = — 0, 
respectively. The other root P = of the saddle point 
equation ( |5q ) is discarded because the prefactor vanishes 
there. In the neighborhood of the saddle points P s , the 
exponent behaves like 



P 2 +c 2 



P 2 + c 2 



(P 2 +2c 2 ) 1 /2 (p2 +2c 2)l/2 

3c 4 



(P 2 + 2c 2 ) 5 / 2 



times ±2ikr. We note 



(P 2 + 2c 2 ) 



P=±iV3c 



P=±j\/3c+e 



He 



P=±iV3c+e 
(P - Ps)' 



-o) 
-o). 



Let us try the right members of the two twins, P s + 
iiv^c + O. Then, 



P + c 6fcr 
± 2ikr— — tttt = t4cA;t ± (P 

(P 2 + 2c 2 )V2 C ^ 



Ps+) 2 , 



(57) 



which requires that the path of integration be deformed 
to run along the imaginary axis and the quadrant of a 
large circle to come back to the real axis. But, on the 
circle, the exponential behaves like 



exp 



±2ikr 



P 2 



(P 2 + 2c 2 ) 1 /2 



exp [T2krRsm(f>] (58) 



(P := Pe i0 ) requiring that < cj> < ir/2 for t > 0. We 
remark that \f(u)p){f(uip) + 1} | behaves on the circle as 
e -ii | cos 20|^ n0 £ decreasing as R — > oo when (j> = ±7r/4, 
while helping suppress (pq) when cf> — 0,ir. 

Thus, we take the path of integration passing the sad- 
dle (§7|), {iy + : y runs from to ± R)} U {Re^ : 
(j) runs from ± 7r/2 to 0)}, with R — > oo, obtaining the 
same results for the two terms ± in (pa): 



(jo(x,t)) 



(c) _ 



.7/2 



2(2^)6 V 25 V3/3c 2 
smn — - — 



x / P a (k) {I 



3/2 



,ik-x— 4ckt j3 



d J k, 



(60) 



which is to be added on (|52 

Note that this result implies a kind of slow oscillatory 
diffusion as illustrated by 

7T\3/2 f , „, 47T 5 / 2 fOO -(ict-ix)k 



(fcr) 



8tt 5 / 2 sin0(]iK|,t) 



{x 2 -l-(4ct) 2 } 



174' ( 61 ) 



where the sin of 



0(1*1.*) := ^ tan ~ 1 ^| 

is responsible for the oscillation mentioned. 

Thus, we have seen that an irreversible behaviour can 
arise from the Hamiltonian dynamics with the initial con- 
dition having the temperature inhomogeneity represnted 
by the density matrix (|2^) and with the model Hamilto- 
nian (|22|) which is time-reversal invariant. 



V. DISCUSSION 

In this paper, we have first pointed out that the thermo 
field dynamics reveals a nonlocal structure hidden behind 
the local-equilibrium density matrix with temperature in- 
homogeneity in the standard statistical mechanics. 

We have then shown with a simple model that, in a 
system of interacting particles, diffusive behaviour can 
arise from the Hamiltonian dynamics with an initial con- 
dition having temperature varying slowly in space even 
when the Hamiltonian is time-reversal invariant. The dif- 
fusive behaviour is seen when we look at the system in 
macroscopic scales both in time and space. The macro- 
scopic space scale is incorporated in our calculation by 
expansion in powers of the wave numbers, and the macro- 
scopic time scale by applying the saddle point method to 
the integration over angular frequency. 

It has to be remarked, however, that the diffusive be- 
haviour (p2|) we have obtained has in x-space the form, 



F W = - J\cV 2 (^) 3/2 /(cpJ{/(c P J + l}e" 4 ^. 



(59) 



We remark that the double poles P = ±i\/2c due to 
f(u>p){f(ujp) + 1} have the residues killed by the expo- 
nential. 

In (p9|), A appears only in e - 4cfe (*+ lA / 3 ) anc j ma y De 
neglected because of the small factor k. Then, 



r 3/2 e -m(x-,) 2 /(«)j( ?/ ) (| 3 yi 



(62) 



multiplied by a damping factor e~ 2mc if written fully 
with the particle mass m, the strength of the interaction 
c and H restored. 

The damping just mentioned and the slow oscillatory 
diffusion we saw in (61) are interesting new features. The 
latter may be a sound wave associated with the diffusion, 
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which could have been smeared out to become invisible 
in macroscopic treatments, while the former could prob- 
ably be an artifact due to our model Hamiltonian not 
conserving the particle number. 

We notice further that the diffusion constant in ([32]) is 
independent of the strength c of the interaction against 
our expectation. This is a result of using expansions of 
different quantities in powers of the wave number k and 
that, as is seen in (45) for example, many of the expansion 
coefficients have the interaction strength c in their de- 
nominators, limiting the validity of our results, say (62), 
to k < c or \x — x'\ > 1/c. Therefore, we cannot let c — > 
without losing the space region of validity of our results. 
If c was 0, in fact, we would lose the diffusive behaviour 
because 9 P should vanish and consequently F' s,± ^ = in 
© and ©. 

After seeing all these features, we stress once again that 
we have shown by a simple example that the Hamiltonian 
dynamics with an initial condition having temperature 
inhomogeneity can lead to an irreversible behaviour even 
when the Hamiltonian is time-reversal invariant. This is 
a fact we have established mathematically, though our 
model and some of its consequences are not quite physi- 
cal. 

The irreversibility so derived stands against Zubarev's 
assertion (|] quoted at the end of Sec. Ill, but not entirely 
against his arguments because he argues simply that the 
results from Hamiltonian dynamics with a time-reversible 
Hamiltonian should be time-reversal symmetric. Our re- 
sult is in fact symmetric under time reversal t — > — t as 
remarked earlier. It is of course not symetric under time 
reversal with respect to any instant of time other than 
t = 0. The instant t = is special because the initial 
condition, which itself is time-symmetric, is given at this 
instant of time. 



[7] 

[8] 

[9] 
[10] 

[11] 



Amsterdam (1982). 

N.N. Bogoliubov, Journal of Physics (USSR), 11, 23 - 32 
(1947). 

D.N. Zubarev, Nonequihbrium Statistical Thermodynam- 
ics, TV. by P.J. Shepherd and ed. by P. Gray and P.J. 
Shepherd, Constants Bureau, New York-London (1974). 
Sec. 20.5. 

K. Huang, and C.N. Yang, Phys. Rev. 105, 767 - 775 
(1957). 

K. Huang, C.N. Yang, and J.M. Luttinger, Phys. Rev. 
105, 776 - 784 (1957). 

C.C. Bernido, M.V. Carpio-Bernido, K.Nakamura and 
K.Watanabe ed., Proc. 1998 Jagna Workshop, to be pub- 
lished, Gordon and Breach, Tokyo. 



This work was presented at the International Work- 
shop on Thermal Field Theories and Their Applications 
in Regensburg, Germany, Aug. 10 - 14, 1998. 



[1] M. Ichiyanagi, Phys. Rept, 262, 227-310 (1995). 
[2] H. Nakano, Int'l. Jour, of Mod. Phys., B7, 2397-2467 
(1993). 

[3] R. Kubo, Jour. Phys. Soc.Jpn, 12, 570-586 (1957). 

[4] H. Ezawa, in Progress in Quantum Field Theory, ed. by 

H.Ezawa, S.Kamefuchi, North-Holland (1986). pp. 305 - 

324. 

[5] H. Ezawa, K. Watanabe, and K. Nakamura, in Thermal 
Field Theories, ed. by H. Ezawa, T. Arimitsu, and Y. 
Hashimoto, North Holland, Amsterdam (1991). pp. 79 - 
94. 

[6] H. Umezawa, H. Matsumoto, and M. Tachiki. Thermo 
Field Dynamics and Condensed States. North Holland, 



8 



